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Abstract 

In this paper we solve the hydrodynamical equations of optically thin, steady state accretion disks around Kerr 
black holes. Here, fully general relativistic equations are used. We use a new method to calculate the shear tensor 
in the LNRF (Locally Non-Rotating Frame), BLF (Boyer-Lindquist Frame) and FRF (Fluid Rest Frame). We show 
that two components of shear tensor in the FRF are nonzero (in previous works only one nonzero component was 
assumed). We can use these tensors in usual transonic solutions and usual causal viscosity, but we derive solutions 
analytically by some simplifications. Then we can calculate the four velocity and density in all frames such as the 
LNRF, BLF and FRF. 



KeyWords: black hole accretion disks, relativistic disks, ac- 
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1. Introduction 

Accretion disks are important in several astrophysical sys- 
tems. They can be found around Young Stellar Objects(YSO), 
around compact stellar objects in our galaxy and around several 
super-massive black holes in Active Galactic Nuclei(AGN). 

In super massive accretion disks the mass of black hole is 
(10 5 — 1O 9 )M0. To study such disks we use the general rel- 
ativity with relativistic hydrodynamics in Kerr metric back- 
ground geometry. In relativistic Navier-Stokes fluid we have 
stress-energy tensor which is related to viscosity and is the 
cause of redistribution of energy and momentum in fluid. This 
tensor is defined by the four velocity and metric. In the previ- 
ous studies, the relativistic and stationary solutions of standard 
black hole disks were solved. Lasota (1994) was the first who 
wrote down slim-disk equations which include relativistic ef- 
fects. He also assumed that only i-<f) component of the stress 
tensor is nonzero. He used a special form for this component, 
which was followed by Abramowicz et al(1996 and 1997). 
Chakrabarti (1996) derived the transonic solutions of thick and 
thin disks for a weak viscosity. He assumed similar form for 
viscosity as Lasota (1994). Then Manmoto (2000) derived the 
global two temperatures structure of advection-dominated ac- 
cretion flows (ADAFs) numerically by using full relativistic 
hydrodynamical equations including the energy equations for 
the ions and electrons. 

Papaloizou & Szuskiewicz (1994) introduced a phenomeno- 
logical and non-relativistic equation for the evolution of vis- 
cous stress tensor (causal viscosity) which has been used by 
many authors. Gammie & Popham (1998) and Popham & 
Gammie (1998) solved ADAFs with relativistic causal viscos- 
ity. They used the Boyer-Lindquist coordinates. 



Takahashi (2007b) solved the equations of relativistic disks 
in the Kerr-Schild coordinates by using the relativistic causal 
viscosity. In the papers of Gammie & Popham and Takahashi 
they assumed that in the FRF, only the r — <f> component of 
shear viscosity is nonzero, then they used the transformation 
tensors to derive the components of shear stress tensor in the 
Boyer-Lindquist or Kerr-Schild frames. 

In present study, we concentrate on the stationary axisym- 
metric accretion flow in the equatorial plane. We use a new 
method to calculate the shear tensor and azimuthal velocity 
of fluids in the locally non-rotating frame (LNRF) by using 
Keplerian angular velocity. We derive two kinds of shear ten- 
sors in LNRF and BLF; in the first one the direction of fluid 
rotation is the same as that of black hole (f2 + ) and in the sec- 
ond one the direction of fluid rotation is in opposite to that of 
black hole(f2~). We calculate the components of the shear ten- 
sor for two kinds of fluids in LNRF, BLF and FRF, these cal- 
culations show that in FRF there are two nonzero components 
(r — t and r — <j> components). But in previous papers, the only 
nonzero component in FRF was r — 4> component. The r — t 
component results from the relativistic calculations of the shear 
tensor which changes some components of the four velocity. 
Then, by using these shear tensor components we calculate the 
four velocity in LNRF, BLF and density in all frames . 

This paper's agenda is as follows. We introduce the metric 
and reference frame in §2. In §3 basic equations are given. In 
§4 the shear tensor is calculated in FRF, LNRF and BLF. We 
derive four velocity in LNRF and BLF in §5 and the influence 
of two important parameters on density and four velocity can 
be seen in this section. Also, the influence of the r — t compo- 
nent of shear tensor can be seen in this section. Summery and 
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conclusion are given in §6. 

2. Metric, Reference Frame 

2.1. Back Ground Metric 

For back ground geometry, we use the Boyer-Lindquist co- 
ordinates of the rotating black hole space time. In the Boyer- 
Lindquist coordinates, the Kerr metric is: 

ds 2 =g afj dx a dx f) = -a 2 dt 2 + -y ij {dx i + P i dt)(dx j 

+/3 j dt), (1) 

where i,j = r,6 and <p. Nonzero components of the lapse func- 
tion a, the shift vector j3 % and the spatial matrix jij are given 
in the geometric units as: 



a = 



"X' 



0* 



7r 



A' 



1ee= S, 



700 : 



Asin 2 9 



(2) 



Here, we use geometric mass m = GM/c 2 , S = r 2 + a 2 cos 2 8, 
A = r 2 -2Mr + a 2 and A = £A + 2mr(r 2 +a 2 ). The position 
of outer and inner horizons, r±, are calculated by inserting 
A = to get r± = m ± (to 2 — a 2 ) 1 / 2 . The angular velocity 
of the frame dragging due to the black hole rotation is w = 
~9t<f>/g4>4> — 2 " ar , where M is the black hole mass, G is the 
gravitational constant and c is the speed of light and the angular 
momentum of the black hole, J, is described as: 



a = Jc/GM 2 



(3) 



where — 1 < a < 1. 

Similar to Gammie & Popham ( 1 998), we set G = M = c = 1 
for basic scalings. The nonzero components of metric, g^, and 
its inverse, g^ are calculated in Appendix 1. 

2.2. Reference Frame 

In our study, we use three reference frames. The first 
one is the Boyer-Lindquist frame (BLF) based on the Boyer- 
Lindquist coordinates describing the metric, in which, our cal- 
culations are done. The second one is the locally non-rotating 
reference frame (LNRF) which is formed by observers with 
a future-directed unit vector orthogonal to t = constant. By 
using the Boyer-Lindquist coordinates, the LNRF observer is 
moving with the angular velocity of frame dragging (u). The 
third frame is fluid rest frame (FRF), an orthonormal tetrad ba- 
sis carried by observers moving along the fluid. 

The physical quantities measured in LNRF are described by 
using the hat such as u* 1 and up and in FRF by using paren- 



theses such as u^) and 



The transformation matrixes of 



FRF, LNRF and BLF are given in Appendix 2 (Bardeen 1970; 
Bardeen, Press & Teukolsky 1972; Frolov & Novikov 1998). 

3. Basic Equations 

The basic equations for the relativistic hydrodynamics are 
the baryon-mass conservation (pu^).^ = and the energy mo- 
mentum conservation; T£ v = 0, where p is the rest-mass den- 
sity and T^ v is the energy-momentum tensor. Basic dynamical 
equations except the baryon mass conservation are calculated 



from the energy-momentum tensor, T^ v . We use the energy- 
momentum tensor written as: 



= pnu^u v +pg^ + iT + q^u v + gV, 



(4) 



where p is the pressure, n = (p + u + p) / p is the relativistic 
enthalpy, u is the internal energy, t^ v is the viscous stress- 
energy tensor and is the heat-flux four vector. The rela- 
tivistic Navier-Stokes shear stress, t^ u , is written as (Misner, 
Thorne & Wheeler 1973): 



t^ 



-2\<j^ - (Oh 



(5) 



where A is the coefficient of dynamical viscosity, ( is the co- 
efficient of bulk viscosity, W" = g^ v + u^u u is the projection 
tensor, = u?^ is the expansion of the fluid world line, and 
a^ u is the shear tensor of the fluid which is calculated as: 



1 



1 



(6) 



where = u^u 1 is the four acceleration. 

We study a stationary, axisymmetric and equatorially sym- 
metric global accretion flow in the equatorial plane, i.e., we 
assume ug = 0. We also assume that the effects of the bulk 
viscosity and heat-flux four vector is negligible. In the follow- 
ing sections, we derive the basic equations using the vertical 
averaging procedures around the equatorial plane which were 
derived in, e.g., Gammie & Popham (1998). 

3.1. Mass conservation 

The equation for the baryon mass conservation is written as: 

1 



(pu") ; „ = -^(V^gpun,^ = o, 



(7) 



where u M is the four velocity and \J~^g = r 2 . By averaging the 
physical quantities around the equatorial plane and assuming 
constant M (the mass-accretion rate ), we can write equation 
(7) as follows: 

(4irH e r 2 pu r ), r = => -ATtH g r 2 pu r = M, (8) 

where Hg is half-thickness of the accretion disk in 9 direction. 
If we normalize the rest-mass density , p, by setting M = 1, for 
calculating global structure of accretion flow, we have: 



-4irH e pu r r = 1, 
also, by differentiating equation (9) we have: 

d\np 



dlnHg 2 d\nu r 
dr r dr 



+ ■ 



dr 



0. 



(9) 



(10) 



3.2. Killing Vector 

Two specific killing vectors of Kerr metric are = (1,0,0,0) 
and = (0,0,0,1), which are used to derive the disk equations. 
At first, angular momentum conservation can be derived by £^ 

{Tffi);* = 0=>>(TZ). v = 0. (11) 
By vertically averaging equation (1 1) we have: 

4r(r 2 7T), r = Mnl - 4irH e r% = Mj, (12) 
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where Mj is the total inward flux of angular momentum. 
This equation is similar to the angular momentum equation of 
Gammie & Popham (1998). 

Another equation which can be obtained by using is: 



(W=0=>(3T)i, = 0. 



(13) 



Similarly, by vertically averaging and inserting constant E, we 
have: 



AirH e r 2 {{p + p + u)u t u r + t r t ) = E. 



(14) 



This equation expresses the constancy of mass-energy flux E 
in terms of the radius; E is the actual rate of change of the 
black hole mass. If the fluid is cold and slow at large radiuses, 
then E w M (Gammie & Popham 1998). 

3.3. Energy Equation 

The equation of local energy conservation is obtained from 
UfjT.^ = 0, as follows: 

r d{p + u) 



dr 



(p+u+p)e + ( P +u+ p)u v u^% 



+t^u„ = 0, 
where is defined as: 

e = «;'„. 

And from equation (7) we have: 



(15) 



(16) 



dp 



{pu v ), v = => p, v u v + pu» v = -j-u r + pu» v = 0. (17) 

Therefore the energy equation can be written as: 

r du u+pdp 



-) = (p + u + p)a"u M + 



v ■ 



(18) 



p dr ' 

In this paper we do not derive the temperature, pressure and 
internal energy so that, we will not use this equation. If we 
want to derive these variables we must use a state equation and 
use the relation of shear tensor components. 

4. Shear Tensor 

Lasota (1994) used the r — <fi component of shear stress ten- 
sor as follows: 



t T 4, = ~VP 



A 3/2 A l/2 7 3 dn 



(19) 

(20) 



r 5 dr 
7 ^(l-(^) 2 )-^, 

where is the azimuthal component of velocity in the LNRF 
which is introduced by Manmoto (2000): 

i A 



r 2 A l/2' 

where f2 is defined as: 

f2 = Q — lu. 



(21) 



(22) 



Abramowicz et al (1996 and 1997), Manmoto (2000) and 
others also used equation (19) for stress tensor. In 1994, 
Papaloizou & Szuskiewicz introduced a non-relativistic causal 
viscosity which was used in relativistic form by Peitz & Apple 



(1997a,b), Gammie & Popham (1998)and Popham & Gammie 
(1998) in Kerr metric and Takahashi (2007b) in Kerr-Schild 
metric. They assumed that in the FRF all components vanish, 
except t r( f, = t<f, r . The t r< p — S can be introduced in relativistic 
form as: 

PS _ S-Sp 

Dt T r ' 



(23) 



where D/Dt = u^Q^ and So is the equilibrium value of the 
stress tensor and t t is the relaxation time scale. Therefore, in 
steady state we have : 



r dS S — Sq 
dr T r 



(24) 



We use a new method to derive the shear tensor approximately. 
We calculate from equation (21) in LNRF assuming u r = 
(the azimuthal velocity is much greater than accretion velocity 
or radial velocity, except very close to the inner edge, and for 
fluids with a small viscosity). In order to simplify, we assume 
= Slfc, where Keplerian angular velocity is defined as: 

" ' (25) 



= ±- 

k r 3 / 2 ±aM 1 / 2 ' 



therefore, for azimuthal velocity we have: 

2ar s 



v r = 



A ±1 

7m\ 



r 2 A l/2V r 3/2 ±a A 

First we calculate the u* = 7 = 



■)■ 



(26) 



i-(t>*) 2 ' 

We can calculate four velocity in LNRF (Appendix 3). The 
four velocity for ft ~ (direction of fluid rotation is opposite to 
the black hole rotation) is: 



ry A(r 2 — a) 



l (r - 2 — a) 2 A — A 2 — Aa 2 r 2 (r'2 — a) 2 — 4arA(r'2 — a) 



,0,0, 



-(r 3 + ra 2 + 2ar^) 



yV 4 (r - 2 — a) 2 A — A 2 — 4a 2 r 2 (r 2 — a)^4arA(r2 — a) 

(27) 

Calculations show that in LNRF = and 6 = 0. The 
shear tensor can be calculated by inserting the four velocity 
of LNRF in equation (46). The nonzero components of u7 

a p,( a (i = U A;7 M ^) anc l f° r ^ m LNRF are: 

ur„ = — -L (+24r 6 a 3 + 12ar 8 - 36r 5 a 3 - 18ar 7 

* ;r 2y/AB- 

+ 12a 5 r 4 + 6a 5 r 3 + + 3a 6 r^ + 2a 6 r^ - 36r^" a 2 

+5r^a 2 + lOr^a 2 + 7r^"a 4 + 4r^a 4 + 4r5a 4 ), 



!!,.= 



(+2r~a 6 + 9ar^ + \2r^a A - 18ar^ 



2y/AB~ 

-3a 5 ri - 2a 5 ri +r 9 + 18r 5 a 2 + 4r 7 a 2 - 20r 6 a 2 
+3r 5 a 4 -4r 4 a 4 -2r 3 a 4 ), (28) 



a- = U? M r = 0, 

t t-r 



a- =u. V + u, ;U = O, 

r r-t r;<p 
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a7=u7 u r = 0, (29) 5. Deriving four velocity 

I Assuming special values for A (for example in Takahashi 

g fr = o H = 4 r^ B - (+ 24r n ' + l2ar -36ra - 18or (2007b), A = 1.5, 1.7, 2.1, 2.2, 2.3, 2.4 and 2.5), we can solve 

19 7 b ii 2 equations (10), (12), (14) and (18) to derive four velocity, den- 

+12a r +6ar + r 2 + 3ar 2 +2ar 2 — 36r 2 a sity, pressure and etc, numerically. In numerical solutions we 

+5r^a 2 + lOr+a 2 + 7r~ a 4 + Ar^a 4 + 4r 2 a 4 ), must insert values of physical variables on boundary conditions 

1 n 1B 9 3 9 specially on horizon. Some of these boundary conditions are 

<J f < j> = = 4y / A_g- (~^^ r 2 a ^ ar 2 ^ r2 a + r nonphysical, therefore we want to derive analytical solutions 

13 7 5 without any boundary conditions. We use some non-relativistic 

-18ar* - 3a 5 r 2 - 2a 5 r% + 18r 5 a 2 + 4rV - 20r 6 a 2 relations to calculate four velocity and density in LNRF and 

+3r 5 a 4 - 4r 4 a 4 - 2r 3 a 4 ), (30) BLF analytically. At first, the t r(f> relation of stress tensor in 

n q 9 s9 Takahashi (2007a) is used in a relativistic disk,this relation is: 

where B = (r 7 — 2r 2 a — r 4 a 2 — 3r b — 6r s a + r b a 2 - 

7 5 3 (]£t 

2r 2 a 3 — 4a 3 r^)3. Also for following Appendix 3, the = -vpr 2 — . (34) 

four velocity is: ^ r 

Similar to Abramowicz et al.(1996) we assume r\ = 1, therefore 
from equation (5) by zero bulk viscosity and A = prjv, we have: 



rv / A(^ 2 + a) 

r 4 (rl + a) 2 A - A 2 - 4a 2 r 2 (ri + a) 2 + AarA{A + a) ~ T dr 7 

therefore, f2 is 



1 ' " <rt> = -zr-r, (35) 



2 — 2ar 2 , ,->+ f 2 



)• ^=l^o%dr. (36) 



r 4 (r 2 +a) 2 A - A 2 - 4a 2 r 2 (rt + a) 2 + 4arA(r 2 + a) 

For calculation the 51 ± we use equations 30 and 33 with a sim- 
^ ' plification in -gq- and -g-r as: 



Therefore, the nonzero components of u^o and for S1+ ^ 

fi-,a A = 0): ^ 

1 7 o o 5 . 3 

1 / * R 51 ^ 8 . nn S 3 . 7 I O^o^.o 1 /l„^o -o 



are(similar to ft", a A = 0): "5T = ( r? + 2r ~ a ~ - ^ ~ 6r3fl2 + r5fl2 

B + 



u+ = — = (-24r 6 a 3 -12ar 8 + 36r 5 a 3 + 18ar 7 +2r 2 a 3 +4a 3 r 2 

* r 2VAB+" 1 3 



(1 + (3r 6 - 2ar~ - r 5 a 2 + r 4 a 2 + 6r 3 a 2 



-12aV 4 -6aV+r- +3aV 2 +2a°r5 -36r~a 2 2r 7 

+5r^a 2 + 10r^a 2 + 7r^a 4 + 4r5 a 4 + 4r5a 4 ), 5 2 n 2 3 ,35. 15 6 11 4 2 

7 ' -r 5 a -2r 2 a-4a 3 r 2 ) + — — (3r b + 2r 2 a + r 4 a 2 

_i_ —1 . _ 11 o _ 15 „ 2 ■< „ „ 12 8r 
uT = — = (— 2r 2 er - 9ar 2 - 12r 2 a + 18ar 2 , „ ». „ 7 , , 5 n9n 

<^ r 2VAB+ +6r 3 a 2 -r 5 a 2 + 2r5 a 3 + 4a 3 r3) 2 ), (37) 

+3a 5 r 3 + 2aM + r 9 + 18r 5 a 2 + 4r 7 a 2 - 20r 6 a 2 
-2r 3 a 4 ), 

1 -2ri a 3 - 4a 3 r3) 



' ' 7 ~ — A ~ R 3„2 I J>2 



+3r 5 a 4 - 4rV - 2r 3 a 4 ), (32) — = (r 7 - 2r 2 a - r 4 a 2 - 3r 6 - 6rV + /■ a 



B~ 

7 - 3 i„S.'-" 



r, =a ■ ■■ (-24r 6 a 3 -12ar 8 + 36r 5 a 3 + 18ar 7 , 3 



r+ - ^+ - 



^ « 4VAB+ V " ~ -7=(1 + A(3r 6 + 2ar^ -r 5 a 2 + r 4 a 2 + 6r V 

,/^2i v 2r 7 v 



-12a 5 r 4 - 6o 5 r 3 + + 3a 6 r 2 + 2a 6 r 5 - 36r~a 2 



+5r^a 2 + lOr^a 2 + 7r^a 4 + 4r§a 4 + 4r?a 4 ), -r 5 a 2 + 2r 2 a 3 + 4a 3 rt) + — (3r 6 - 2r~a + r 4 a 2 



15 
87 1 

0-+, =<rt = — — (-2r^a 3 -9ar^ - 12r^a 3 + r 9 +6r 3 a 2 - r 5 a 2 - 2r 2 "a 3 - 4a 3 ri) 2 ). (38) 

r<t> <t>r 4 V / A_B+ 

+18arf + 3a 5 ri + 2a'V t + 18r 5 a 2 + 4rV - 20r 6 a 2 InfluenCeS ° f M& fPP™^ in shea ^ ^ components 

4 4 4 4 are seen in figure 1 for ir and a = .9. Solid curves have no 

+3r a — 4r a — 2r a ), (33) simplifying and dotted curves are with this simplifying. 

where B+ = (r 7 + 2r^a - r 4 a 2 - 3r 6 - 6r 3 a 2 + r 5 a 2 + 

2r 2 a 3 + 4a 3 r 2 The shear tensor in BLF can be derived c ... , r> 2 rr 1 u u- <->nm \ t u 1 

- I . ,. , ~ , - . . Setting / = S2r z (Takahashi 2007a), the angular momentum can 

from cr a g — e'ie R or l y in which ef; and e« are the transformation , , , . , c t . , , , . , 

, p . 01 p >z, rrr, " r be calculated. From equations (36) and (12), p can be derived 

tensors(Appendix 2). There are two nonzero components of 

shear tensor in BLF, er tr and <T r( ^, which are given in Appendix 

4. Also using transformation matrices, two nonzero compo- p _ jjli j (39) 

nents of shear tensor in FRF can be calculated including &(t)(r) 9,-KHgr 2 va^ r 

and a {r)w which are given in Appendix 4 as well. ^ y . g (Takahashi 2007a) . 




Fig. 1. Influence of this simplifying in shear tensor components for Q + and a - 
curve are with this simplifying 



.9. Solid curves are without that simplifying and dotted 



v = aa 2 jn k , 



(40) 



which is the usual a prescription for viscosity(Shakura & 
Sunyaev 1973). In some papers He is calculated by study- 
ing the vertical structure such as Riffert & Herold (1995), 
Abramowicz et al.(1997) and Takahashi(2007b). In these pa- 
pers He was calculated by vertical averaging, then introduced 
by other variable such as I, u t , p, p and etc. If we use each of 
those relations of He, we must solve all equations numerically. 
But we want to derive an analytical solution therefore similar 
to Takahashi(2007a) we use a simple form of He as: 



He 



5 a s 
2 fifc 



(41) 



definitely, p is the same in LNRF and BLF and is shown in 
figure 2 for and (£! + ) by assuming a s = .1. 



Also u r can be calculated from equation (9) and Uf can be 
found from equation (14). Assuming E ~ M = 1 we have: 



1 



-E--^ 

pu r 



(42) 



47T/W * 

We can calculate four velocity in BLF using the transformation 
equations in Appendix 2, as: 



(43) 
(44) 




Now, we may see {\u t \ , |u t -|), (\u r \, \u f \) and Qu^ 
figure 3. 



(45) 
)in 



Four velocity and |0| in BLF(|u e = 0|) are shown in figure 4 
setting a s = .1. 



In equation (40), < a < 1 is viscosity coefficient in a pre- 
scription which influences on redistribution of energy and mo- 
mentum. The influence of these parameters are shown in figure 
5 with a, = .1 



The four velocity and density can be derived for various total 
inward flux of angular momentum (for example j = 1,2 and 3) 
which are shown in figure 6. 



5.1. Influence ofr — t component on the four velocity 

In pervious papers such as Abramowicz et al. (1996 and 
1997), Gammie & Popham (1998), Manmoto (2000) and 
Takahashi (2007 a,b), only the r — <p component of shear tensor 
was assumed to be nonzero in FRF, but in section 4 we showed 
that the r — t component is also nonzero. The r — t compo- 
nent of shear tensor in LNRF results from covariant derivative 
of tt£. In general relativity the gravitating field causes time di- 
lation. Due to this time dilation, the coordinate time(t) and 
proper time(r) are not the same(dt > dr), therefor u l in LNRF 
can be derived as: 



dt 
dr" 



(46) 



tt* = |u t -| is different in each radius and also tangent of Uf cre- 
ates the r-t component of shear tensor. To see the influence of 
this component, we derive four velocity with both r — t com- 
ponent of shear tensor and without it. If we set of = in the 
pervious section the figures are similar to that of the previous 
papers. From equations (39) and (42) it is clear that r — t com- 
ponent affects m in all frames because without this component, 
the second term of u t in equation (42) vanishes. If we calculate 
w/ 1 = g^ v u v , the r — t component also affects in BLF. The 
influences ofr — t component of shear tensor are shown in fig- 
ure 7, solid curves are with this component and dotted curves 
are without it. 
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Fig. 3. Four velocity in LNRF and BLF. Left: fi~ , right: fi+ (a = .01, j = 3, a = 



— .9 and a. 



= .1, solid: in BLF and dotted: in LNRF) 
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Fig. 4. Four velocity and |fi| in the BLR Left: Q , right: fi + (a = .01, j = 3, solid: a = .9, dotted: a = .4, dash-dotted: a — 
long-dashed: a = — .4 and dash-spaced: a = — .9) 
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Fig. 6. Influence of j in density and the four velocity of 51 + in BLF(q = .01, a = .9 and a s = .1. Solid: j = 1, dotted: j = 2 and 
dash-spaced: j = 3) 



No. ] 
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Fig. 7. Influence of the r — t component of shear tensor on |it t -j(in LNRF) and |w t |, |tt*|, u^(in BLF), |fi| and |w r |for solid with the 
r — t component and dotted without it(a = .01, a = .9, o s = .1 and j = 3). 



10 



Moeen Moghaddas et al. 



[Vol. 



Obviously, r — t component of shear tensor have more influ- 
ence on ut and u*. Therefore, for deriving ut and u* we can 
not vanish the r — t component. 

6. SUMMERY AND CONCLUSION 

In causal viscosity method, only r — (f> component of shear 
tensor is assumed to be nonzero in FRF. We do not use this 
method for viscosity because our calculations show that there 
are two nonzero components of shear tensor in FRF. In our 
method, we use azimuthal velocity of LNRF and Keplerian an- 
gular velocity, then we calculate all components of shear ten- 
sor in LNRF for two kinds of fluids (rotation in the direction of 
black hole and rotation in the opposite direction of black 
hole(fi~)). Using transformation tensor we can calculate all 
components of shear tensor in BLF and FRF. 

Solid curves of figure 1 show unphysical treatments in shear 
tensor components close to the inner edge. This unphysical 
treatment in a^ f may be concerned to assuming u r « u^, 
therefore if we find a suitable v r this treatment may be re- 
solved. We may have a suitable v r if we put w* = 7 = 

y Vt^-t^) 2 > ^ 1 ~ ( ^ )2 > (l/)2 ' Accordin 8 t0 ec i ua - 

tion (21), the most important origin of unphysical treatments in 
ct.^ is the 1/A term (1/A has a singularity in the inner edge). 
Therefor by using Taylor expansion the unphysical treatment 
of uif. and were solved as it can be seen in dotted curves of 
figure 1. 

We solved the hydrodynamical equations in LNRF with as- 
suming 77 = 1, we derived the density and four velocity in 
BLF and LNRF. The density which is calculated analytically 
in LNRF is the same in all frames such as BLF. But com- 
parison of four velocity in LNRF and BLF (figure 3) shows 
that u r and u t have a sm all difference in tw o frames. In 
equation (44), ^/S/A s=s \Jr 2 j (r 2 — 2r + a 2 ) is near 1, spe- 
cially in larger radiuses, therefore u r is similar in two frames. 
Following equation (43), for u t we have two parts, the first 
part is (y/AE/A w y(r — 2)/r)u t - and the second part is 
(— 2arj\jATj w — 2a/r 2 )u^. Because value of is greater 
than Uf (dotted curves in figure 3), the first term has more influ- 
ence on Ut than the second term (the second term is important 
just in inner radiuses). 

and are different in two frames, because in equation 

(45), y/ A/H s=s r. The previous equation is logical because the 
LNRF rotates with the angular velocity of the frame dragging 
(the frame dragging velocity can not be seen in LNRF but it can 
be seen in BLF), therefore is greater than in all radiuses. 

In section 5, four velocity was calculated analytically in 
LNRF and BLF. The sign of £1 introduces two types of flu- 
ids, + and Q~, but it can be seen that for example, in figures 
2 or 4, in a = .9 f2 + is similar to 0~ in a = — .9 (and in a = .4 
Q+ is similar to Q~ in a = —.4 and so on) therefore we use 
in figures 5, 6 and 7. 

The effects of a coefficient are shown in figure 5. It shows 
that if a increases, then the shear viscosity grows, the radial 
four velocity will increase and then density decreases. 

The r — t component of the shear tensor in LNRF which re- 
sults from the relativistic calculations(equation 46)is equal to 



covariant derivative of u t ~. If we ignore the r — t component 
of shear tensor (the dotted curves in figure7), the figures are 
similar to the figures of the previous works such as Gammie & 
Popham (1998) and Takahashi (2007b). But, with r - 1 compo- 
nent of shear tensor (the solid curves in figure 7) a few number 
of four velocities change in all frames, specially u*, u*, and 
v r will change more than the others. The influences of the r — t 
component on all physical variables can be derived if we solve 
all the equations of the disk with a state equation numerically. 

When we have time dilation u l has to be greater than l(dt > 
dr u l — > 1). According to figure 4, the effect of time 
dilation is much greater in the inner edge than the other places. 
If we ignore the r — t component of shear tensor or put \u^\ = 
l(dotted curve in the first panel of first column of figure 7), it 
is equivalent to ignoring time dilation. Equations (42) and (44) 
show that the r — t component of shear tensor have no influence 
on u r and u r . Near the black hole u l increases, which causes a 
decreases in v r (v r — \). 

Energy can be calculated from T u — prji/u 1 +pg u , where 
the first term of energy is related to u t . Figure 2 shows that p is 
greater in inner edge than the other places and also is the same 
with the r — t component of shear tensor and without it. The 
Ut which uses the r — t component of shear tensor is greater 
than the u t without it. Therefore, the first term of energy with 
r — t component of shear tensor is greater than the this term of 
energy without r — t component. 

Far from the black hole the general relativistic influences are 
vanished, therefore in outer edge, the influences of r — t com- 
ponent of shear tensor are too low (as it is seen in outer radiuses 
of figure 7). 

If we want to calculate the temperature, pressure, internal 
energy, cooling and heating or radiation, we must use a suitable 
state equation with equations (9), (12), (14) and (18)(energy 
equation), then, to solve the equations numerically suitable 
boundary conditions are needed. 
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Appendix 1. Metric Components 

Nonzero components of Kerr metric are given as: 

„2 i a . Q<t> 
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and inverse components of metric are: 
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Appendix 2. Transformation Between BLF, LNRF and 
FRF 

Components of connecting between BLF( Boyer- 
Lindquist Frame) and LNRF(Locally non-rotating Frame) are 
calculated as: 
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The transformation between LNRF and FRF(Fluid Rest Frame) 
are as follows: 
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Where we use u — —u^ = au , then the Lorentz factor 7 is 
calculated as: 



7^a-« 2 r 



2 = cm 



(v 2 = v l v t = vj + vj+v 2 -).(Al) 



In LNRF, three velocity components are calculated as: 

v l = — (i = r,6,4>). 
After some calculation we have (Q — ^): 



(A8) 
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Appendix 3. Deriving Four Velocity in LNRF 



(A3) From equation (26) we have: 
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First we derive four velocity for flk + = l/(r* + a) as 

(All) 



^ 2 A 2 +4a 2 r 2 (r3 +a) 2 - 4arA(r^ + a) 
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In LNRF 7 = u* = y j-^p- and in our study we put tig = 0. 
If we suppose >> u? (as it was in pervious papers such as 
(A4) Popham & Gammie 1998 and Takahashi 2007a,b), then v 2 w 
(u^) 2 and we have: 

7 = u* = — ti£ = ry/A(r i + a) x 

(r 4 (r^ +a) 2 A- A 2 -Aa 2 r 2 {ri + a) 2 + AarA(ri +a))~i, 

(A12) 

also we know that = u 4 ju l therefore, 
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?/ = u*V = (r 3 + ra 2 - 2ar'i ) x 

(r 4 (ri +a) 2 A-A 2 -4a 2 r 2 (ri + a) 2 + 4arA(ri + a))~5. 

(A13) 

For 07/ ) = — 1/ (r 5 — a) we have: 
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Appendix 4. Nonzero components of shear tensor in BLF 
and FRF 



Therefore, in FRF, two components of shear tensor are non- 



zero. 



Nonzero components of shear tensor in BLF can be cal- 
culated by a a p = e^epO-ftv where and e v are given in 
Appendix B. Therefore, these components for Q + are: 
1 



3r 2 a 6 + 8r 8 a 3 + 36r~a 2 + 4r^a 6 



4B+A-5 

24r^a 4 - 4r 4 a 5 - 40rV - 8r 5 a 5 + 2r 10 a + 6r 6 a 5 



-18r-a 2 - 4r^a 4 + 36r 6 o 3 ) 



+ - 



1 



-(+3a 8 r^ +r¥ + 42r n a - 12a 7 r 6 + 2r ^a 8 



6r 5 a 7 



4B+A2 

+6r* a 2 + 12r^ a 4 - 12r 12 a - 36a 5 r 8 + 54rV 
+12a 5 r 6 - 36r 10 a - 40a 4 r^ - 72r 8 a 3 - 36r 10 a 3 - 2r* 
-56r^a 2 - 8r^a 4 + lOr^a 6 + 72r^a 2 + 102rV), 

4. 4. VA 



^ _-:3« '/- • 2a '/- • 3a r I)/- a 

+18r 5 a 2 - 12r^a 3 + 18r^a - 2r^a 3 - 2r 3 a 4 - 20r 6 a 2 
-4r 4 a 4 + 4r 7 a 2 + r 9 ). (A15) 

Non-zero components of shear tensor in FRF can be calculated 

by o-(a)(/3) = e \ a ) e \^) u ^ where e^ a) and efo are transforma- 
tion matrixes of Appendix B. 
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Also, we have: 



